Bifurcation Phenomenon in a Spin Relaxation 
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Spin relaxation in a strong-coupling regime (with respect to the spin system) is investigated in 
detail based on the spin-boson model in a stochastic limit. We find a bifurcation phenomenon 
in temperature dependence of relaxation constants, which is never observed in the weak-coupling 
regime. We also discuss inequalities among the relaxation constants in our model and show the 
well-known relation 2IY > I\, for example, for a wider parameter region than before. 



Nowadays, we are gaining technologies for control- 
ling quantum states, and there are increasing necessities 
for better understandings on dissipation or relaxation in 
nanoscale systems. Many authors have been working on 
this kind of issue, and we already have some knowledge 
about it, especially in the field of quantum optics, where 
the weak-coupling limit, e.g., the van Hove limit [Q, §, 
is a very good approximation. There are, however, many 
other situations we have to handle. 

In this Letter, we focus our attention on a dissipative 
dynamics derived through a different scaling limit from 
that of van Hove's; the limit A — > with the scaled time 
t = X 2 t fixed is applied to a total Hamiltonian of the 
form §, | 



fltot = A 2 ff s + XV + H E 



(1) 



where H$ is the Hamiltonian of a dissipative system of 
interest, Hb that of an environment with infinite degrees 
of freedom, and XV the interaction Hamiltonian. In con- 
trast to the van Hove limit ||, §], which is applied to a 
Hamiltonian of the form H tot = Hs + XV + Hb, the limit 
for the Hamiltonian (Q) is regarded as a strong- coupling 
limit in the sense that the interaction XV is much larger 
than the energy scale of the system X 2 H$- 

The systems with the scaling ([|) have been discussed 
by several people ||, and their importance is recog- 
nized already. The present authors have recently con- 
structed p] a framework of the stochastic limit for the 
scaling modifying the original stochastic limit which 
is proposed by Accardi et al. || |(| on the bassis of the van 
Hove limit. Applying the framework for the scaling (|l|) 
to the spin-boson model Q, which is simple but is fun- 
damental and has many applications, we here investigate 
the spin relaxation under the scaling ([!]) in detail and 
report interesting features in the relaxation dynamics, 
which are never found in the weak-coupling regime. 

We consider the spin-boson model 



H S = -cr z + —<r x , 



V = iha, 



H R = 



(2a) 



du> 



(gMou-g'Wat) (2b) 



under the scaling ([!]) . Here a x , er 2 are the Pauli matrices, 
and the spin system Hg has two eigenstates 



Hs\±) 



±-Ml \±) 



m = Ve 2 + a 2 



(3) 



The environment Hb is the boson system, where a w (aj,) 
is the annihilation (creation) operator of the boson of 
mode uj. Two systems interact with each other through 
the linear coupling interaction V . In the following anal- 
ysis, the form factor g(u>) is assumed to be such that 
2ir\g(ui)\ 2 ~ rjuj for to ~ (Ohmic model). We assume 
further that, before time t = 0, the spin and boson sys- 
tems are uncorrelated and the boson system is in the 
thermal equilibrium state at temperature T, i.e., the ini- 
tial state of the total system is given by 



Ptot(0) = As ® p B , 



PB = |e 



-H B /k B T 



(4) 



where ks is the Boltzmann constant and Z the normal- 
ization constant. 

Under these assumptions, it is possible to derive a mas- 
ter equation for the density operator of the spin system, 
ps (t) = trs ptot (t) . In the stochastic limit for the scal- 
ing (Q), the master equation reads j| 



^Ps(r) = - 1 -[Hs,ps{t)] -\\<t„ [<t,,Ps(t)]\, (5) 



where 



7 



g _ 2rjk B T 



(6) 



The equivalent Bloch equation for the spin system op- 
erators D± = |±)(T|, Do = |+X+| - |-}(-| is given, 
through the relations trs[/9s(r)Z)±] = trs[psD±(r)], 
^s[ps{t)D ] = tr s [psD (t)}, by 



dr 



D{t) = -AD(t), 



(7a) 



where 
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(7b) 



and e = e/hilo, A = A/HUq. Based on the master equa- 
tion ((H) or the Bloch equation (Q), we investigate a spin 
relaxation in the scaling ([!]). 

The dynamics of the spin system is characterized by 
the eigenvalues of the matrix A in the Bloch equation (j^) . 
There are three eigenvalues, which are analytically ob- 
tained by Cardano's formula for the solutions to third 
order polynomial equation. In some cases, two of them 
are complex conjugate to each other, and the remaining 
one is real. The real eigenvalue gives an exponential de- 
cay, and the complex ones damped oscillations (Fig. [I]). 
We call the former "longitudinal relaxation" and the lat- 
ter "transverse relaxation." The damping rates T l for 
the longitudinal relaxation and Tt for the transverse one 
are defined, respectively, by 



= (real eigenvalue), 

= (real part of the complex eigenvalues). 



(8) 



Our definitions of the "longitudinal" and "transverse" 
relaxations seem formal, and the physical meaning of 
them is not so clear as is in the weak-coupling regime. In- 
deed, in the latter case, the matrix A in the correspond- 
ing Bloch equation is always diagonal A — diag(IY — 

iVl Ri t t + in R ,r L ) with r L = AV, r T = AV/2 

[7 9 = 27r| 5 (Oo)| 2 coth(^o/2fcsT 1 ) and n R is a (renor- 
malized) frequency] ||, and the longitudinal relaxation 
is nothing but the energy relaxation. In the Bloch equa- 
tion (0), Tl stands just for a damping rate of an ob- 
servable (not necessarily energy) which decays exponen- 
tially without oscillation, and it depends on tempera- 
ture (see below). Furthermore, there are cases where 
all three eigenvalues are real giving three "longitudinal" - 
relaxation constants (i — 1,2,3). It might be inap- 
propriate, therefore, to introduce such notion as "longi- 
tudinal" or "transverse" in this regime. One should no- 
tice, however, that a quite nontrivial phenomenon, i.e., 
bifurcation in the relaxation dynamics, can be observed 
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FIG. 1: (a) Longitudinal and (b) transverse relaxations 



in terms of F^ and Tt as is seen in the following. We 
are interested, not in the matter of definition, but in the 
fact that the phenomenon peculiar to nonlinear (and/or 
nonequilibrium) systems can be seen in this system. 

Now, let us investigate the relaxation constants T t and 
Tt in detail. First of all, it should be mentioned that pos- 
itivity of the relaxation constants is guaranteed since the 
master equation (||) is of the Lindblad form || |9) . For an 
investigation of the matrix A in the Bloch equation (Q), 
the following proposition will be useful: 



Proposition 1 Let 

M 3jR = {A : 3 X 3 matrix| tr A, det A, tr adj A e 



(9) 



where adj A is the cofactor matrix of A (adjugate of A) 
and M is the set of real numbers. The necessary and 
sufficient condition for the real parts of the eigenvalues 
Xi (i — 1, 2, 3) of A G M% :R to be positive semi-definite 



ReA,>0 (i = l,2,3) 



(10) 



is that 



trA > 0, det A > 0, tradjA > 0, and f(trA) > 0, 



(11) 



where /(A) = det (A - A). 



(Its proof is quite elementary and straightforward jHl.) 
The matrix A in the Bloch equation (]?]) belongs to M^ tR 
and satisfies the condition (11). The real parts of the 
eigenvalues of A, i.e., the relaxation constants, are there- 
fore positive semi-definite by Prop, [j] pi. 

Next, let us discuss the parameter dependence of the 
relaxation constants Tt and Tt- Since we have derived 
the master equation (||) or the Bloch equation (]?]) from 
the microscopic point of view, i.e., from the basic Hamil- 
tonian (|l|) with (Q), the relaxation constants Tt and Tt 
are given in terms of the parameters in the Hamiltonian, 
e, A, •q 1 and the temperature T . Furthermore, since the 
relaxation constants are obtained analytically, their de- 
pendence on the parameters is clear. In Fig. |[ temper- 
ature dependences of the relaxation constants Tl and 
Tt are presented for various values of e/KIq. The first 
point to notice is that they exhibit quite different be- 
haviors depending on the parameter e/MIq. This is not 
the case in the weak-coupling regime, where the relation 
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FIG. 2: Temperature dependences of the longitudinal- 
relaxation constant Tl (solid lines) and the transverse one 
T T (dashed lines) for (a) e/hQ, Q = 0.0, (b) 0.2, (c) 0.4, and 
(d) 0.6. 
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FIG. 3: Characteristic functions f (A) = det(A — A) in the 
temperature regions I— III in Fig. 0(b). 



2Tt = Tl always holds. In particular, there are tem- 
perature regions where no transverse-relaxation constant 
exists [Figs. ||(a) and (b)]. 

Furthermore, observe an interesting and nontrivial 
phenomenon there. See Fig. ||(b), for example. A 
transverse-relaxation constant bifurcates into two longi- 
tudinal ones at some temperature (at C from lower tem- 
perature and at A from higher temperature). Or con- 
versely, two longitudinal ones merge into one transverse 
one (at C from higher temperature and at A from lower 
temperature) . 

The mechanism is understood by tracing the zero 
points of the characteristic function /(A) = det(A — A), 
i.e., the real eigenvalues of the matrix A. The forms of 
/(A) in the temperature regions I— III in Fig. ||(b) are 
sketched in Fig. |L The real eigenvalue in region I disap- 
pears in region III and changes into two complex eigen- 
values when the temperature becomes higher. On the 
other hand, a new zero point comes out in region II, i.e., 
two complex eigenvalues in region I change into a real 
one. These are typical bifurcation mechanisms, which 
are often seen in nonlinear (and/or nonequilibrium) sys- 
tems [ |l2| |. This is a bifurcation phenomenon seen in a 
spin relaxation process. 




FIG. 4: Temperature dependence of the longitudinal direction 
I corresponding to Fig. 0(b). In region II, only two directions 
which connect continuously with that in region I or III are 
displayed. 



An experimental implication of this phenomenon is as 
follows. Let us focus on a spin component 



0~L 



I ■ (T = l x a x + La y + l z a 



(12) 



which decays exponentially without oscillation: o~l(j) = 
o~lz~ VlT ■ The vector I is normalized to unity, and we 
call it "longitudinal direction." Temperature dependence 
of the longitudinal direction I is shown in Fig. ||[ corre- 
sponding to those of the relaxation constants in Fig. ||(b). 
In the limit T — > 0, I is in the "direction of the Hamilto- 
nian" (er^ cx Hs), and in the limit T — > oo, it is in the 
"direction of the interaction" (ox oc V). If one traces I 
from T = to higher temperature, it disappears at some 
critical temperature [A in Fig. ^ or |2|(b)] and is found 
in a different direction B. If one traces it from higher 
temperature in reverse, it jumps at a different tempera- 
ture from C to D. The bifurcation phenomenon will be 
observed experimentally as a discontinuous temperature 
dependence of the longitudinal direction. 

Let us finally discuss the ratio of the longitudinal- 
relaxation constant to the transverse one Tt- Tem- 
perature dependences of it are plotted in Fig. ||, corre- 
sponding to Fig. 0. They imply that the inequality 



2T T > r L 



(13) 



holds at any temperature T and for any e/hTlo, while the 
equality 2Tt = Tl holds irrespectively of the parameters 
in the weak-coupling regime ||. Actually, the inequal- 
ity ([[3]) is proved on the basis of the master equation (j|) 
or the Bloch equation (0). For that purpose, the follow- 
ing proposition is available: 



Proposition 2 If A s satisfies the condition (p 

then the necessary and sufficient condition for the set of 
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FIG. 5: Temperature dependences of Tl/Tt, corresponding 
to Fig. ^. Note that there are regions where no transverse 
relaxation constant exists. 



inequalities 



Re Ai + Re A 2 > ReA 3 , 
Re A 2 + Re A 3 > ReAi, 
Re A 3 + Re Ai > Re A 2 



(14) 



(H is an Hermitian operator in 2-dim. Hilbert space and 
{cij} is a positive semi- definite c- number matrix), and 
showed that the inequalities (|l^) and ([l6]) hold if the 
superoperators TL and V commute: 



[H,V] = 0. 



(18) 



In our master equation (g) , which is also of the Lindblad 
form, they do not necessarily commute. (Consider, for 
example, the case of T > and e ^ 0.) This means that 
the inequalities (O) and ( |l6| ) hold more generally f p"3[ . 
Notice that Eq. (|l3| ) is a well-known inequality in actual 
experiments 



14 



It is true, of course, that the model Hamiltonian (|l|) 
with (||) is to be refined for a description of real ex- 
periments, but we expect the interesting and nontriv- 
ial features shown here to be universal ones which re- 
flect some important characteristics related to nonlinear 
and/or nonequilibrium systems. 
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One of the authors (K. I.) is grateful to Luigi Accardi for 
kind hospitality. He is supported by an overseas research 
fellowship of Japan Science and Technology Corporation. 



to hold is 



/(trA/2) >0. 



(15) 



(Its proof is also elementary and straightforward H.) 
Since the matrix A in the Bloch equation ^ satisfies the 
condition (15) as well as (ll]), the relaxation constants 
satisfy the inequalities (|l4) by Prop. [| When two of 
the three eigenvalues of A are complex and the rest is 
real, the inequalities ( |l4| ) are nothing but the inequal- 
ity (p"3|). Furthermore, the above proposition has shown 
that when there exist three longitudinal-relaxation con- 



stants (i = 1, 2, 3), they satisfy the inequalities (14) 
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The inequalities (|l^) and ( |l6| ) are also discussed by 
Gorini et al. |)| : They investigated the Markovian master 
equation of the general form, i.e., the master equation of 
the Lindblad form M fif, for two- level systems, 



-p(t) = (H + V)p(t), 



(17a) 



Hp=--[H,p], 



Vp 
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<lin{[<Ji,p<Jj\ + [0-iP,(Tj]} 

(17b) 
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